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Abstract 

We consider the branching and annihilating random walk A — ► 2A 
and 2A — * with reaction rates a and A, respectively, and hopping 
rate D, and study the phase diagram in the (X/D,a/D) plane. Ac- 
cording to standard mean-field theory, this system is in an active state 
for all a/D > 0, and perturbative renormalization suggests that this 
mean-field result is valid for d > 2; however, nonperturbative renor- 
malization predicts that for all d there is a phase transition line to an 
absorbing state in the (X/D, cr/D) plane. We show here that a simple 
single-site approximation reproduces with minimal effort the nonper- 
turbative phase diagram both qualitatively and quantitatively for all 
dimensions d > 2. We expect the approach to be useful for other 
reaction-diffusion processes involving absorbing state transitions. 
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1 Introduction 



Branching and annihilating random walks (BARW pQ) have been the focus 
of much attention [2J El HI El El E] j as they are among the simplest models 
of nonequilibrium critical phenomena observed in physics and other sciences 
(for reviews see, e.g., [3EKTU]). They are generic reaction-diffusion processes 
in which particles of some species A move stochastically on an arbitrary d- 
dimensional lattice and are subject to the creation and annihilation reactions 

A ^> (m+l)A and kA — ^> 0. The nonequilibrium phase transitions of these 
models are known to belong to two distinct universality classes, depending on 
the parity of m and k: the "Parity Conserving" class (for m and k even) and 
the "Directed Percolation" class (for m odd) p2, E] • We are interested here 
in BARW belonging to the Directed Percolation class, which are generally 
denoted "odd-BARW" . 

The simplest odd-BARW, which in this work we will call for short the 
"OB A model" (odd branching and annihilating walks), is defined by 

' A^2A 

OBA model I 2A ^ (1.1) 
k ^ U 

where a and A are on-site creation and annihilation rates, respectively, and 
D is a hopping rate between adjacent sites. If the lattice has coordination 
number c, this means that a particle will leave its lattice site at a rate cD. 
Since it captures the essential critical properties of its entire class [Sj, our 
focus will be specifically on the OBA model. Although it is well-established 
that the OBA model is in the Directed Percolation universality class [5], it 
has appeared much harder to determine its phase diagram, and this is the 
subject of this paper. 



In standard mean-field (MF) theory [H E] the OBA model is described 
by the rate equation 

DAp + ap- Ap 2 , (1.2) 



dt 

where p(r, t) > is the particle density. For all branching ratios o > 
Eq. (11.21) has two spatially uniform stationary solutions, the "absorbing state" 
Po = and the "active state" p* = ct/A. For a > the active state is stable 
and is reached exponentially fast in time from any initial state with p > 0. 
For a = the OBA model coincides with the pure pair annihilation process 
(PA) and corresponds to a critical point in parameter space at which the 
absorbing state is reached according to the power law decay p(r,t) ~ t~ x . 
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On a lattice in finite dimension d the question of the stationary states 
of fll.ip is much more difficult to answer. For o = 0, i.e. for the PA, 
perturbative renormalization group analysis has shown [TTJ [5] that below the 
critical dimension d c = 2 the decay to the absorbing state slows down due 
to depletion zones (anticorrelations) in the spatial density distribution and 
follows the power law p(r,t) ~ t~ d l 2 . 

Early simulations of the OBA model (3J, H] showed that for d = 1 and 
d = 2 absorbing states exist even in an interval of branching ratios a > 
0, in contradistinction to MF theory. This indicates, therefore, that low- 
dimensional fluctuations qualitatively change the "phase diagram" of this 
system in the (A, a) plane. 

In a seminal paper, Cardy and Tauber [5] formulated a field theory for 
general BARW, which they analyzed by perturbative renormalization group 
techniques. For the OBA model they concluded that in dimension d < 2 & 
minimum branching ratio a c is needed in order for the system to be able to 
sustain an active state. For small A the critical value a c behaves as 

ff - B (2iLf- A ^°- d=2 - e - (L3) 

The analysis by Cardy and Tauber is valid for small A and appears to break 
down for d > 2 |5J . However, since A becomes irrelevant above two dimensions 
and since from the PA analysis one can assume that fluctuations are also 
small for d > 2 in the OBA model when a is small, Cardy and Tauber argue 
that MF theory should be restored for d > 2, that is, the system should be 
active for all a > [5]. 

This picture was modified by an analysis due to Canet et al. P [7J, who 
employed nonperturbative renormalization group (NPRG) techniques. They 
demonstrated that in all finite dimensions d > 2 there is a threshold A c such 
that 

(i) for A < A c the MF result of an active state for all o > remains 
true, but 

(ii) for A > A c there is an active state only when a > cr c (A). 
Hence according to the work of Ref. [H [7J standard MF theory fails dramat- 
ically for this system. 

Following this failure, the OBA model has been re-investigated [13] by 
means of an alternative MF type approximation, namely the "cluster MF 
method", also called "generalized MF method" and originally proposed for 
non-equilibrium systems in [12]. This approach consists in considering the 
master equation for blocks of iV sites and truncating the hierarchy of prob- 
ability distributions, so that it can be solved numerically. From cluster MF 
calculations of the OBA model Odor [13j has confirmed for N > 2 the ex- 
istence of a finite threshold A c > above which an inactive phase exists. 
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The purpose of this paper is to set up a simple analytically tractable 
approximation that correctly predicts the existence of an absorbing phase 
in all finite dimensions in the appropriate domain of the parameter space. 
We propose a new approach, to be called single-site approximation, which 
allows diffusion steps to take place only to empty sites. The underlying idea 
is that in dimensions d > 2 where the intersection of two directed random 
walks becomes unlikely, the destruction mechanism that drives the system 
to an absorbing phase is dominated by "on-site" annihilation of a particle 
with its own offspring on the same site, and that we may neglect the annihi- 
lation of particles that meet due to random diffusion. Hence the single-site 
approximation has the nature of a tree diagram approximation in coordinate 
space. It has the virtue that it leads to what is essentially a single-site prob- 
lem which allows for analytic results to be easily obtained. For the OBA 
model these results closely reproduce the NPRG phase diagram, including 
its dimensional dependence, for all dimensions d > 2. 

We believe that in future studies this method may find application to 
other reaction-diffusion processes having absorbing states and that it may 
provide essential information about their phase diagrams ahead of any more 
sophisticated work on them. 

In section [2] we define the single-site approximation for the OBA model, 
which for short we will refer to as the "SS-OBA model." We obtain its 
solution in section 12.21 We show in section 12.31 that the model tends to an 
absorbing state in a specified region of the (A/-D, o~/D) plane whose shape we 
discuss in section 12.41 Our results are qualitatively and quantitatively fully 
consistent with the NPRG results. In section [3] we add various comments to 
the discussion. Section 0] is a brief conclusion. 

2 Single-site approximation 
2.1 Definition 

We define the SS-OBA model as follows. The stochastic motion of its parti- 
cles is governed by the rules: 

(i) Each particle is subject to the on-site creation reaction A — > 2A at 
a rate a. 

(ii) Each pair of particles on the same site is subject to the annihilation 
reaction 2 A — > at a rate A. 

These two rules are therefore the same as in the original OBA model. 

(iii) Each particle may hop away from its site at a rate cD and always 
arrives on an empty lattice site (of some abstract lattice that need not be 
specified - say for instance to one of the next nearest neighbours if all neigh- 
bouring sites are occupied). 
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This rule differs from the corresponding one in the OBA model; it means 
that in the SS-OBA model the notion of lattice structure is lost. 

For D = the two models are identical; we therefore expect the SS-OBA 
model to be a good approximation of the OBA model in the small diffusion 
regime. 

2.2 Solution 

To see that the SS-OBA model is exactly soluble, it suffices to note that 
no particle ever enters an occupied site from the outside and that therefore 
each active site has a dynamics independent of the others; a site's occupation 
number evolves only due to on-site creation and annihilation transitions and 
to departures. The solution therefore decomposes into the analysis of the 
time evolution on a single site given its initial condition at some time to, and 
the analysis of the coupling between sites due to the diffusion mechanism. 
These two questions are studied in subsections 12.2.11 and 12.2.21 respectively. 

2.2.1 Single-site problem 

Let P(n, t) be the probability that a specific site contains exactly n particles 
at time t. Then this probability satisfies the master equation 

a(n - l)P(n - 1, t) - anP(n, t) 
+ ±A(n + l)(n + 2)P(n + 2,t)- \\n{n - l)P(n, t) 
+cD(n+l)P(n+l,t) - cDnP(n,t) (2.1) 

for n = 0,1,2,... and with the convention that P(—l,t) = 0. Here the 
lattice coordination number c is the only parameter reminiscent of the original 
lattice. We introduce the scaled variables 

r = At, a = (r\-\ D = cDX' 1 (2.2) 

and set P(n,t) = p(n, r). Defining the vector p(r) = (p(0, r),p(l, r), . . .) we 
may write ( 12. ip as 

^ = Mp, (2.3) 



dt 



P(n,t) 
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where M is the tetradiagonal matrix 



M 



/ D 1 

-D-a 2D 3 

a -2D -25-1 3D 6 

2a -3D -3a -3 AD 

3a -4L>-4a-6 . 

V ••• ••• / 

(2.4) 

Eqs. f l2.3|) - fl2.4p constitute a problem with two parameters, D and a. Except 
when a = D = we expect Eq. (12. 3 p to have only a single stationary solution, 
viz. p st (n) = <5 n .o- The reason is that as n gets large, the annihilation rate 
dominates the creation by one order in n, which prevents "escape" of the 
site occupation number to n = oo; the occupation number, therefore, can 
get caught only in n = 0, even though for large a it may be in a long-lived 
"metastable" state. When a = D = the particle number can change only 
by pair annihilation, which is a parity conserving process. There will then 
be two independent stationary states, viz. p s q(u) = 5 nj0 and pf(n) = S nj \. 

We will write G nno (r — r ), where r — r > and n, n = 0, 1, 2, . . ., for the 
solution of Eq. f)2.3p with initial condition p(n, r ) = 5 nno . It is not possible 
in the general case to write this solution in an explicit closed form, but we 
will suppose that all its essential properties can be determined. In particular 
we assume that, except when a = D = 0, the function G nn i{j) tends to zero 
exponentially at some time scale l//j,i(a,D), 

G nn iT) ~e" WT , r^oo, (a, D) ^ (0, 0), n,n' > 0, (2.5) 

with ni(a,D) > 0. When a = D = we have /Xi(0,0) = 0, which signals 
the degeneracy of the stationary state. This completes the discussion of the 
single site problem. 



2.2.2 Coupling between sites 

Let at time r = the initial state be such that there are S n (0) sites with 
occupation number n, where n = 1,2.... The total initial particle number 
iV(0) is then given by 

oo 

iV(0) = 5>S n (0). (2.6) 

n=l 

We are interested in the average number (N(t)) of particles at some arbitrary 
instant of time r > 0; here (...) denotes an average with respect to the initial 
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distribution of the S n (0) and the stochastic time evolution. We will proceed 
by first calculating the averages (S n (r)). 

There are two types of sites, those that are occupied initially, and those 
that get occupied only later during the time evolution due to diffusion steps. 
We denote the contribution of these two types of sites by superscripts (0) 
and (1), respectively, so that 

(S n (r)} = {S n (T))V> + {S n (T))V>. (2.7) 

Upon considering the time evolution of the initially occupied sites we find 

oo 

(Sn(r)) {0) = ]TGW(t)M0). (2.8) 

n'=l 

Throughout the time interval < r' < r new occupied sites are created due 
to diffusion steps at a rate D(N(t')). Hence 

POO 

(S n (r))W = D dr'G nl (T-r')(N(T')), (2.9) 



where the upper boundary of the integral has been sent to oo exploiting that 
Gnn'it) vanishes for t < 0. Summing Eqs. (12.81) and ( 12. 9ft yields 

oo 

(S n (r)) = X)GW(T)Sn'(0) 

n'=l 

poo 

+ D / dT'G nl {T-T'){N{r')). (2.10) 

J 

When multiplying this equation by n and summing on n we obtain 

oo oo 

(iV(r)> = ^^^(r)5 rf (0) 

n'=l n=l 

/•OO 00 

+ D / dr'J2^G nl (r-r')(N(r')}, (2.11) 
Jo i 



71=1 



which is a closed equation for (N(r)). 

For convenience let us now restrict the initial states to those that have 
only singly occupied sites; that is, we take S n (0) = N(0)5 n i. This clearly 
does not affect the long time behaviour of the system so it does not restrict 
the generality of our argument. Furthermore we define 

oo 

H(T-T')=J2 nG nl(T-T'). (2.12) 

n=l 
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When substituting the previous definitions in Eq. (12.111) we find 

poo 

(N(t))=N(Q)H{t) + D / dr'H(T - t')(N(t')). (2.13) 



In terms of the Laplace transforms 

/•oo poo 

N(s) = / dr e^ T (iV(r)), H(s) = / dr e~ ST H(r), (2.14) 
Jo Jo 

it becomes N(s) = N(0)H(s) + DN(s)H(s), whence the solution 

N(s) = N(0) — , (2.15) 

l-DH(s) 

which may be inverse Laplace transformed to (N(t)). This completes the 
solution of the average total particle number in the SS-OBA model. 

2.3 Existence of an absorbing phase 

Because of Eqs. (I2.12p and (12.51) the decay of H (r) will be characterized by 
the same /ii as that of G nn i{r). Let us suppose that H{t) = e _AtlT , so that 
H(s) = l/(/ii + s). Then Eq. (12151) implies that N(s) = iV(0)/(//i - D + s), 
whence 

(N(t)) = N(0)e~ { ^~ D)T . (2.16) 
Therefore the condition for (N(t)) to tend to zero is 

Hi(a,D) - D > 0. (2.17) 

The SS-OBA model has an absorbing phase in the region of parameter space 
where Eq. (I27T71) holds. 

To analyze this equation we consider first the special case D = 0. In this 
case the SS-OBA model is described by a single-site master equation identical 
to Eq. (12. ip but with D = 0. It reaches the absorbing state exponentially 
an asymptotic rate /^i(5", 0) which for < a < oo satisfies 

//i (or, 0) > 0. (2.18) 

We note that /Xi(0,0) = as emphasized in section 12.2.11 and that we must 
furthermore have lim^oo yUi(o", 0) = 0, since the decay of the metastable state 
becomes infinitely slow in that limit. We invoke now continuity of fii(a,D) 
in its second argument and conclude that for all values of a in (0, oo) there 
exists a positive threshold value T> c (a) such that Eq. (I2.17P is satisfied for all 
< D < V c (a), i.e. the stationary state of the SS-OBA model is absorbing 



S 



for all < D < T> c (a). Reverting to the original parameters this means that 
for all ratios < cr/A < oo there exists a V c (a/X) such that for 



A, 
D 



> 



V c {a/X) 



(2.19) 



the stationary state is absorbing. Written this way, inequality ( 12.191) allows 
for comparison with the phase diagrams of Ref. [7j which are plotted in the 
plane of abscissa X/D and ordinate a/D and are displayed in Fig.[TJ First of 



d=1 6=2 d=3 d=4 d=6 



a 




Figure 1: Phase diagram of the OBA model (jl.ip ensuing from the nonperturba- 
tive renormalization group (NPRG) in dimensions d = 1, . . . , 6, from Canet et al. 
[7|. For each dimension, the active phase lies on the left and the absorbing phase 
on the right of the transition line. For d > 2 this line comes into the X/D axis at 
some threshold value X c /D with a finite slope onprg ~ 2.3. The spacing between 
the phase transition lines for two successive dimensions is Anprg ~ 2.2. 



all Eq. (12.191) implies the existence for large enough A of an absorbing region 
in the phase diagram, which is in full agreement with the NPRG results of 
Fig.[U Furthermore, as the spatial dimension d tends to infinity, so does 
the lattice coordination number c (typically linearly with d) , and Eq. ( 12.191) 
shows that the region of phase space to which our proof applies, recedes to 
infinity as c —>■ oo. This, too, corroborates the NPRG results of Ref. [7] 
which indicate that an absorbing phase exists in all finite dimensions. It 
also matches with standard MF theory, which has effectively d = oo, and no 
absorbing phase in this limit. 

One can take a further step and analyze the shape of the phase transi- 
tion line between the active and absorbing phases obtained in the SS-OBA 
model. As shown in the next section, it fits with the NPRG results even on 
a quantitative level. 
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2.4 Analysis of the phase diagram 

In this section we analyze in greater details the location of the phase tran- 
sition line in the (X/D,a/D) plane defined by Eq. (I2.19p . We show that it 
intersects the X/D axis at some finite value X/D = X c /D, and with a posi- 
tive slope. To do so we determine //i(cr, D) for arbitrary D perturbatively in 
small a. To linear order in a, the result, derived in the Appendix, is 

/n(a, D) — D — a + 0{d 2 ). (2.20) 

Combining this expression with criterion ( 12.17ft we see that for a — > 0, the 
system will tend towards an absorbing state at the condition that 

D<1 or A>A C , (2.21) 

where 

X c /D = c = 2d, (2.22) 

and in which the last equality is for the case of a hypercubic d-dimensional 
lattice. 

The variation of the threshold X c /D with the dimension has also been ob- 
tained [7] within NPRG. Indeed, from Fig. [I] which is for hypercubic lattices, 
this variation appears to be linear, X c /D ~ A NPRG <i with the slope equal to 
Anprg ~ 2.2 (as estimated in [7j). This result is in very close agreement 
with our Eq. ( 12.221) . which yields A = 2 for the same slope. 

The second order correction in a to /ii can also be worked out, as shown in 
the Appendix, and allows for the determination of the slope of the transition 
line. To second order the smallest eigenvalue is given by 

/ ~s ~ D-l 2(2D 2 -3D-3) 2 
/xi(cr, D) = D — = = L a 2 + 0(a 3 ). 2.23 

According to the criterion of Eq. (I2.17P the phase transition line between the 
active and the absorbing phases is defined by the condition yL\— D = 0. Sub- 
stituting (12.231) in this condition, dividing out a, writing a = (a/D)(D/c), 
and using that a, D — 1, and X/D — c are of the same order we find 

a = (ao+3)(o + i)' D _ l) + 0((6 _ m 

D 2D(2D 2 - 3D - 3) V ' 

UX -c) + o((±-c)'), (2.24) 



2 \D J \\D 

which is the equation for the phase transition line near the threshold point 
(A/D,<7/£>) = (c,0). 
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Figure 2: Phase diagram of the OB A model (jl.ip according to the single-site 
approximation of this work, for the dimensions d and d + 1. The circle at infinity 
has been represented by a dashed arc. The phase transition lines come into the 
X/D axis with a slope a = | and for two successive dimensions have a spacing 
A = 2. This phase diagram is very close qualitatively and quantitatively to the 
one of Fig.CQ 

The slope a at the threshold appears to be a = §, independently of c, 
i.e., of the dimension d. The value a — ~ compares favorably with the NPRG 
phase diagram of Fig. [I] where the phase transition lines appear to be merely 
drifting as the dimension grows, with a quasi-constant slope «nprg numeri- 
cally equal to «nprg ~ 2.3 [14|. 

Fig. [2] summarizes the results of this section and shows that our single- 
site approximation captures all the essential features of phase diagrams of 
the OB A model in d > 2. It allows in particular to probe the dependence on 
the dimension, which is beyond the scope of standard MF approach. 

3 Discussion 

In addition to the discussion that has accompanied the above determination 
of the phase diagram, several points deserve some comments. We provide 
them here. 

Order parameter. Our first remark concerns the active state. According 
to Eq. (I2.17P this state is characterized by the inequality \i\ — D < 0. In 
the region of the phase diagram where this inequality holds, Eq. (12.161) shows 
that the total particle number increases exponentially in time. Nevertheless, 
in this regime the average number of particles per occupied site, (n) occ , is 
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well defined in the limit r — > oo. One may consider this quantity as the order 
parameter of the SS-OBA, but it should be emphasized that the usual order 
parameter is instead p = q(n) occ where q is the fraction of occupied sites. 

In the SS-OBA model each diffusion step takes the diffusing particle to an 
empty site. For the discussion of the phase diagram in the preceding sections 
there has been no need to specify whether this is a site that has perhaps 
previously been occupied or whether it is an entirely new site. Therefore the 
fraction 1 — q of empty sites, and hence the usual order parameter p, remain 
undefined in the SS-OBA model. 

Hence a calculation of p within our approach would require further elab- 
oration and/or modification of the model. A way to do this was pointed out 
by Dickman [15] , but goes beyond our more restricted purpose of analytically 
studying the phase diagram. 

Reformulation. The SS-OBA model may be formulated in an equiva- 
lent way which preserves the original lattice structure. This reformulation 
requires that we distinguish between two notions that coincided in the def- 
inition of section (12. lft . viz. between a "lattice site" in its usual sense and 
a "family" of particles: a family consists of a particle having diffused to a 
(now not necessarily empty) lattice site, together with all the offspring it 
has generated on that site and which has never left it. Initially all particles 
on the same site are considered to constitute a family. Consequently, the 
particles on each site may at any time be partitioned into families. When a 
particle performs a diffusion step, it leaves its family and starts a family of 
its own on its arrival site, where other families may or may not already be 
present. 

We may then replace rules (ii) and (iii) of section 12.11 by the following: 
(ii') Each particle can annihilate (at a rate A per pair) only with other 
members of its own family on that site. 

(iii') Each particle performs diffusion steps to neighboring sites with a 
rate D per transition; for coordination number c this means that a particle 
will leave its site at a rate cD. 

This reformulation of the SS-OBA does not change the mathematics and, 
in particular, leads to the same (N(t)) as found in section |2~2"1 It has the 
merit of bringing out clearly that two-particle annihilation in the SS-OBA 
model occurs under more restrictive conditions than in the original OBA 
model. This makes it tempting to believe that the average total particle 
number in the SS-OBA model is an upper bound to the same quantity in 
the OBA model. If true, our calculation would constitute an exact proof 
of the existence of an absorbing state in the OBA model. We have not, 
however, been able to prove this upper bound property and leave it as an 
open problem. 

Below two dimensions. Our final remark concerns what happens when 
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the spatial dimension d is equal to or less than the critical dimension d c = 2 
for pure pair annihilation. For d < d c the renormalisation group (both pertur- 
bative and nonperturbative) predicts that in the OBA model the threshold 
A c vanishes. In the SS-OBA model this would correspond for instance to 
^c(O) — > oo in Eq. (12.191) for d < 2. We have not investigated this point fur- 
ther, since the behaviour of the system for small X/D and o~/D corresponds 
to the large diffusion regime, for which we do not a priori expect the SS-OBA 
model to be a good description of the OBA model. 

4 Conclusion 

We have formulated a new approximation for reaction-diffusion problems 
with an absorbing-state transition. Its characteristic feature is that it for- 
bids annihilation reactions when one or more of the participating particles 
have moved from the site where they were originally created. The approach 
then leads to what is essentially a single-site calculation, which allows the 
consequences to be determined analytically. 

We have applied the approximation to the OBA model A —>■ 2 A and 2A — ► 
0. We have shown that our theory produces qualitatively and quantitatively 
the main properties of its phase diagram in agreement with the predictions 
of nonperturbative renormalization but with far less effort. 
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A Appendix 

We wish to calculate the smallest nonzero eigenvalue, called —fi\(a,D), of 
the matrix M defined by Eq. (12.41) . We will perform this calculation for 
arbitrary D and perturbatively for small a. The symbol 1 will denote the 
identity matrix. For any matrix L with rows and columns labeled by n = 
0,1,2,..., we write for the matrix obtained from it by erasing its rows 
and columns of indices n = 0,1, . . . ,j — 1; and we denote by [L] mm / the matrix 
obtained from it by erasing the row m and column m! . Hence M^ ^ = M and 
M« = [M] 00 . 
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If we modify the matrix M by suppressing its subdiagonal (of elements 
0, a, 2<7, 35", . . .), it becomes upper triangular. This modified matrix has 
eigenvalues — 1/& given by its diagonal elements, i.e., 

v k = kD + ka + \k(k - 1) = uf ] + fccr, k = 0,1,2,.... (A.l) 

Hence its smallest nonzero eigenvalue is V\ = D + a. Restoring now the sub- 
diagonal will change the and we will calculate this change perturbatively 
to second order in a. That is, we will look for a solution \i — n\ of Eq. (1A.5I) 
which has the form 

fjLi = V! + Avx (A.2) 

with 

Ai/i = u[ 1] a + 4 2) ^ 2 + 0(a 3 ). (A.3) 
The nonzero eigenvalues fi of M satisfy 

det (M (1) = 0. (A.4) 

A cofactor expansion of (IA.4I) along the column n = gives 

- (^i - fi) det (M (2) + /il (2) ) - 5 det ([M (1) + /il (1) ]i ) = 0, (A.5) 

whence, after we substitute ( 1A.2I) in ( 1A.5I) . 

. det ([MW + gilWjio) 

The first order correction i/^ follows from linearizing Eq. flA.6j) in cr, which 
amounts to setting //i = = D and cr = in the two determinants of 
the right hand side of (1A.6j) . Both determinants then reduce to a product of 
diagonal elements and identical factors cancel. Using the explicit expression 
of M together with (|A.3D we then get from (1A.6I) the first order coefficient 

v? = -JUL, (A.7) 
D + l 

which leads to the result shown in Eq. (12.201) . 

The second order in a in the expansion ( 1A.3I) determines the slope at 
which the phase transition line comes into the X/D axis in the phase diagram. 
This second order correction requires that one computes the determinant 
ratio in Eq. (IA.6I) to linear order in a. Performing a cofactor expansion of 
both determinants along their column n — 0, one gets 

det ([M (1) + / u 1 l (1) ] 10 ) = 2D det (M (3) + /iil (3) ) 

-2 ( Tdet([[M«+ /Ul l( 1 )] 10 ]io), 

det(M^+^) = -(u 2 -^)det(M^ +/i 1 l (3) ) 

-2 ( rdet([M( 2 )+/i 1 l( 2 )] 10 ). (A.8) 



14 



To obtain the expansions of Eqs. (1A.8I) to first order in a, it suffices that 
we evaluate the two determinantal coefficients of a in the second and fourth 
line above to zeroth order. Furthermore, it turns out that the first order 
contributions of det (M^ 3 ) + /ill®) in Eqs. (1A.8I) cancel out in the ratio 
(IA.6P , so that only the zeroth order of this determinant is actually needed as 
well. The zeroth order of all the determinants is again obtained by setting 
/i! = and a = 0, upon which the matrices become upper triangular and 
the determinants reduce to the products of diagonal elements. Finally, 

u? = 2(2P2 :W (A.9) 
(2D + 3)( J D+1) 3 



By combining (1A.2I) . (1A.3I) . (1A.7I) . and (1A.9j) one obtains the second order 
result ( 12.231) exploited in the main text. 
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